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The report is devoted to the resuhs of the numerical study of the virtual cath- 
ode formation conditions in the relativistic electron beam under the influence of the 
self-magnetic and external axial magnetic fields. The azimuthal instability of the 
relativistic electron beam leading to the formation of the vortex electron structure 
in the system was found out. This instability is determined by the influence of the 
self-magnetic fields of the relativistic electron beam and it leads to the decrease of the 
critical value of the electron beam current (current when the non-stationary virtual 
cathode is formed in the drift space). The typical dependencies of the critical current 
on the external uniform magnetic field value were discovered. The effect of the beam 
thickness on the virtual cathode formation conditions was also analyzed. 
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I. INTRODUCTION 



The analysis of nonlinear microwave oscillations and complex structure formation mecha- 
nisms in spatially extended systems with intensive beams of charged particles in the regimes 
of the virtual cathode (VC) formation attracts great attention of scientific community^"-. 
It is well known^^Aio that the systems with VC are characterized by the complex dynamics 
and can demonstrate a wide range of nonlinear phenomena, including dynamical chaos and 
pattern formationii"— . Microwave generators using electron beams with VC (vircatorsiii^) 
are perspective devices of high-power microwave electronics for the generation of the im- 
pulses of wide-band microwave radiation due to its high output microwave power, a simple 
construction (particularly vircators can operate without external focusing magnetic field), a 
possibility of a simple frequency tuning and regime switching^i^ii^"— . All these circumstances 
increase the fundamental and applied importance of studies of the nonlinear dynamics of 
the electron beams with VC. 

The oscillating VC is known to appear in an electron beam when the beam current ex- 
ceeds a certain critical value, Icr (space charge limiting (SCL) current and the beam 
space charge is strong enough in order to form a potential barrier (VC) which reflects the 
electrons back to the injection plane. The mechanisms of the VC formation have been inves- 
tigated in detail in case of the one-dimensional (ID) electron beam motion (fully magnetized 
beam)^ii^i^"— , with the critical beam current value being analytically defined for this case 
in Ref.i^. 

However, the use of ID-theory to study the vircator systems is inefficient in many cases 
because it ignores a lot of important factors in the VC behavior and does not agree often 
with the experimental results. So, the development and use of 2D and 3D models for the 
analysis of the dynamics of electron beam with VC have attracted the great interest of 
scientific community recently^'2^~— . The important problem in this field is the analysis of 
the VC formation conditions and nonlinear dynamics of electron structures in relativistic 
electron beams (REB) in the presence of finite external magnetic field or even without 
external magnetic field which focuses the REB. In particular, the studies of REBs with 
overcritical currents are necessary for the analysis of contemporary high-power devices with 
VC - relativistic vircators^ and ion acceleration systems^i^^. 

Analyzing REBs, it is necessary to take into account effects being insignificant for weakly 
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relativistic beams, in particular, the influence of the self-magnetic field of the REB that 
effects considerably on the system dynamics in case of ultra-relativistic electron beams. 
Therefore, the 3D fully electromagnetic self-consistent model of REB dynamics is required 
for accurate and correct analysis of the VC formation in this case. 

The present report deals with the 3D numerical electromagnetic study of the VC non- 
linear dynamics of the annular REB in the presence of an external finite uniform axial 
magnetic field. The structure of this paper is the following. Section [III contains the brief 
formalism describing a 3D mathematical model for the nonlinear interaction simulation of 
electron beam with overcritical current and electromagnetic fields. Section IIIII deals with 
the nonlinear dynamics of REB with overcritical current. We analyze the influence of the 
external magnetic field on the VC formation conditions and structure formation in REB 
and discuss the azimuthal REB instability which leads to the formation of a vortex electron 
structure. In Section HVl the influence of the beam thickness on the VC formation conditions 
is analyzed. The conclusions of this paper are summarized in Section |Vl 

II. GENERAL FORMALISM 

The model under study consists of finite-length cylindrical waveguide region (electron 
beam drift chamber) with length L, radius R and grid electrodes at both ends. An axially- 
symmetrical monoenergetic annular electron beam with the current /, electron energy We-, 
radius Rb and thickness d is injected through the left (entrance) electrode. Electrons can 
leave the waveguide region by escaping through the right (exit) grid or by touching the 
side wall of the drift chamber. In the present work the values of geometric parameters 
were chosen following: L = 40 mm, R = 10 mm, R^ = 5 mm, d = 1.5 mm (except for the 
Section IIVI where the influence of the beam thickness on the VC formation conditions is 
analyzed) . 

The external uniform magnetic field with induction Bz = Bq is applied along the longitu- 
dinal axis of waveguide. The electron beam source is supposed to be magnetically unshielded 
in the considered model. This assumption means that the external magnetic field in the drift 
tube is equal to the magnetic field in the electron source region, therefore, the electron beam 
doesn't acquire azimuthal velocity components at the injection plane (in accordance with 
Busch's theorem21). Such magnetic field distribution is typical for many devices of the 
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high-power electronics, particularly for magnetically isolated diodes that forms high-current 
REBs^I. 

Time-dependent fully 3D electromagnetic model of REB dynamics based on solving the 
self-consistent set of Maxwell equations and equations of charged particles motion accompa- 
nied by corresponding initial and boundary conditions (particle-in-cell method)^!^ is used 
in the present paper. The main equations are written as: 



IdU 19E Air. 



^ = Ei + |pi, BJ/7,, ^ = Pi/7i, i = l...A', (2) 

where E and H are the electric and magnetic intensities, p and j are the charge and current 
densities, r, p, 7 = (1 — (f/c)^) v are the radius vector, impulse, relativistic factor and 
velocity of the charged particles, correspondingly. The subscript i denotes the number of 
particle and is the full number of particles used to simulate the charged particles beam. 

The numerical simulation 3D scheme is based on 2.5D model developed in our previous 
work—. The equations of charged particles motion are used for electron beam simulation 
and solved numerically by means of Boris algorithnt^. The longitudinal Vz, radial Vr and 
azimuthal vg velocity components are calculated with the help of this algorithm on each 
time step. 

The electromagnetic fields in the drift chamber of REB are obtained by means of the 
numerical solution of the Maxwell's equations ([1]) in cylindrical geometry on the shifted 
spatio-temporal meshes with constant spatial longitudinal, Az, radial, Ar, and azimuthal, 
AO, steps and time step, At^^'^^. The values of steps of spatio-temporal meshes is picked out 
from the Courant-Friedrichs-Lewy condition^^i^. The every field component is calculated on 
the own mesh (see^i^). The space charge and current densities on the meshes are calculated 
using a bilinear weighing procedure^^. To model the electromagnetic power output we use 
the approach^ based on the filling the section of electrodynamical system {L < z < 1.2L) 
with the conducting medium with the conductance a. 
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III. NONLINEAR DYNAMICS OF RELATIVISTIC ELECTRON BEAM 
WITH OVERCRITICAL CURRENT 



A. Virtual cathode formation conditions 



Let us consider the resuhs of numerical simulation of VC formation features in the annular 
REB in the presence of the external magnetic field. Here we study the infiuence of the 
magnetic field on the critical current of REB (i.e., VC formation conditions) as well as on 
the nonlinear dynamics of the electron structures in the REB. External magnetic field, Bq, 
varies within the range [0, 40] kGs. 




/;ghz 



FIG. 1. Amplitude spectrum of oscillations of self-electric field longitudinal component at the 
VC region for / = 17 kA, Bo = 



VC in the electron beam with overcritical current is characterized by the complex non- 
stationary spatial-temporal oscillations^'^^ that leads to the excitation of vircator's electro- 
dynamic system, and, consequently, to the generation of high-power microwave radiation 
in such system. Fig. [1] demonstrates the numerically obtained typical amplitude spectrum 
of oscillations of self-electric field longitudinal component at the VC region in the consid- 
ered vircator system. The carried out electrodynamic analysis has shown that frequencies 
of spectral components in Fig. [1] are determined by the corresponding eigenmodes of the 
electrodynamic system formed by the drift chamber walls. REB critical current is supposed 
to correspond to such value of the beam current when the spectral components appear at 
the output spectrum of vircator system and electrons start to refiect back to the injection 
plane z = 0. 
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Fig. [2^ shows the dependencies of REB critical current, /cr, on the external magnetic field 
value. Bo, for different beam energies. We- The insert panel in Fig. [2^ demonstrates the 
corresponding critical current values for weakly relativistic electron beam with the initial 
electron energy We = 79keV (see Ref.— i^). Analyzing this Figure, one can see that the 
curves leriBo) for REB have characteristic feature in comparison with the weakly relativistic 
case when the critical current decreases monotonously approaching the minimal constant 
value with the growth of Bq. 
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FIG. 2. (a) Dependencies of the REB critical current on the external magnetic field induction 
value for the following electron beam initial energy values We'- curve 1 corresponds to 480 keV, 
2 - 600 keV, 3 - 850 keV and ^ - 1 MeV. The insert panel demonstrates the dependency of the 
beam critical current on the magnetic field Bq in the weakly relativistic case (beam energy We = 
79keV)^i^. (6) Illustration of the typical REB curve Ict^Bq) with denoted character values of the 
external magnetic field, where is the value of the magnetic field when the critical beam current 
reaches the first minimum; flmax is the magnetic field when the critical beam current reaches the 
local maximum and i?^m magnetic field when the curve IcriBo) saturates 

The typical shape of REB critical current dependency on external magnetic field value 
is shown in Fig. [2]d. As one can see from Fig. [2], the REB critical current increases with 
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the growth of induction, Bq, within the range B < -B^ax- With further growth of magnetic 
field for Bq > -B^^^. the REB critical current decreases and for Bq ~ -Bmm ^^e saturation 
at the constant level is observed. Electron beam appears to be fully magnetized in the case 
of the strong external magnetic field Bq > -B^j„ and moves generally one-dimensionally. 
Bogdankevich and Ruhadze were shown^^ that the analytical expression for the critical 
current value in the case of ID motion of electron beam has the following form: 



where R is the radius of cylindrical drift tube, Rh d are the radius and thickness of annular 
REB, 7o and Vq are the relativistic factor and velocity of REB electrons at injection plane, 
c is the speed of light, and 1] is the specific electronic charge. Table H] demonstrates the 
comparison of REB critical current values obtained analytically using relation with 
the results of numerical simulation in the case of the strongly magnetized beam {Bq = 
40 kG) for the different values of the electron energy We- Analyzing Table [H one can 
see that the numerical results for REB critical currents agree accurately with the values 
obtained analytically, and, therefore, we can conclude that the developed numerical 3D fully 
electromagnetic model gives correct results. Note also, that the numerically obtained critical 
current value exceeds slightly the analytical results. It's a consequence of the 2D dynamics 
of the electron beam occurring in the system even in the case of the strong external magnetic 
field, e.g. the cyclotron rotation of electrons, pulsations of the beam boundary et al. These 
2D effects reduce the charge density in the VC area and, consequently, enlarge weakly the 
critical current value in comparison with ID estimation ([3]) which doesn't take into account 
these effects. 

It should be noted that the effects of the decrease of the space charge limiting current 
in comparison with ID electrostatic case ([3])™ due to the electromagnetic transients in 
the injected currents defined by the inductive voltage, Ldl/dt, and the power loss in the 
drift chamber through the open boundary have been reported earlier— i^^. However, these 
phenomena do not play the significant role for the system under consideration. Actually, 
the inductive transient effect is pronounced only for the short drift chambers {W/L > 1, 
where W = 2i?f, is the beam diameter, and L is the drift chamber length)^, while for the 
considered system W/L = 0.25. The power loss is also negligible in our simulations, so its 
influence on the critical current value is insignificant in the considered case^. 




To = 
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7] d/Rb + 2 \n{R/Rb) ' 
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TABLE I. Critical current values obtained analytically using relation ([3]) third column) and 
numerically (I^"*", fourth column) for the different values of the beam energy {We, first column); 
second column shows corresponding values of vq/c; Bq = 40 kG. 

WekeV, vo/c I«"kA^ 47"^ kA^ 

480 0.8571 4.82 4.93 

600 0.8882 6.45 6.61 

850 0.9270 10.01 10.12 

1000 0.9412 12.11 12.23 



^ Analytical results^^ 
^ Numerical results 



Thereby, the carried out analysis has shown (Fig. [2^) that REB critical current curves 
have local maxima at the certain external magnetic field value Bq = B^^^^ dependent on the 
beam energy Wg. Such behavior is not observed in the weakly relativistic electron beams 
and develops with the electron beam energy growth. 



B. Structure formation in REB with overcritical current 

The observed behavior of the REB critical current IcriBo), (see Fig. [2^) is determined by 
the influence of the self-magnetic fields of REB on the space charge dynamics and electron 
structures formation. Therefore, this effect is more pronounced in the case of the weak 
external magnetic field and high energy of the injected electron beam We > 600 keV, when 
the influence of the self-magnetic flelds of REB is stronger. Let us consider the effect of 
electron structure formation in the presence of the REB self-magnetic fleld in detail. 

Figures |3H5] show the typical phase portraits of the electron beam which are presented by 
the projections of the instantaneous positions (black dots in flgures) of charged particles of 
the beam at the longitudinal {z, r) and transverse planes (r, 6) for the beam currents being 
close to the critical values Ic-^Bq) (/ > Icr{Bo)) and the different characteristic values of the 
external magnetic fleld, B^. Fig. [3] corresponds to the case Bq < B}^^^ when the external 
magnetic fleld value is less than magnetic fleld -B^ax corresponding to the local maximum 
on curve 1 in Fig. |2^, Fig. |4] — to the case of Bq ~ B]^^^ and Fig. [5] — to the case when 
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the magnetic field corresponds to the range of saturation of Jcr (-Bo )-curve on the constant 
level (i.e., Bo > B^J. 
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FIG. 3. Projections of tfie instantaneous positions of the electron beam charged particles at the 
plane {z,r) (left figures) and plane {r,0) for Zs = 6 mm (right figures) at the consecutive time 
points ti and t2 (^2 — = 0.1 ns) for i?o = 3 kGs, / = 7.5 kA; We = 480 keV. Only particles behind 
the projection plane {z > Zg) are shown in Figures. The vertical dashed line with coordinate Zg in 
figures a and c denotes the projection plane for figures h and d) 

The space charge dynamics of the ultra-relativistic electron beam at the VC region differs 
considerably in comparison with the weakly relativistic case for the weak external magnetic 
field < -Bo < -B^^^ (see Fig. |2b). This major difference is determined by the effect 
of the instability of REB found for the wide range of beam currents in the case when 
< -Bo < -B*„„^. The observed instability, the so-called azimuthal instability of REB, leads 
to the axial symmetry loss of the electron beam dynamics. 

This azimuthal instability of initially axially symmetrical beam in cylindrical drift cham- 
ber facilitates the formation of the vortex electron structure in the drift space. This insta- 
bility is caused by the influence of the self- magnetic fields of the REB (see Fig. |3]). Actually, 
there is the intensive transverse beam current in the system as a result of the beam di- 
vergence due to Coulomb's repulsion forces action in the presence of the weak external 
magnetic field which does not restrain these forces. This transverse current results in the 
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FIG. 4. Similar to Fig. [3] but for the following parameters: Bq = 5kGs, / = 9kA; We = 480 keV; 
t2-ti = 0.2 ns 



appearance of the longitudinal self- magnetic field, -Bf, that causes the azimuthal Lorentz 
force action on the electrons moving in the transverse direction?^. Therefore, these electrons 
get the azimuthal velocities and REB starts to rotate around symmetry axis of the drift 
space, r = 0. Centrifugal force acts on rotating electrons and, consequently, the vortex 
electron structure is formed in the beam that results in the strong azimuthal asymmetry of 
the REB. The formed vortex structure rotates in the drift tube as can be clearly seen by 
comparing the configurational portraits in Fig. |3]a and c for two consecutive time moments. 
The physical mechanism of the azimuthal instability development has similar features with 
the well-known convective instability of the electron beam in the finite external magnetic 
field^>^ that arises due to the non-uniform distribution of the beam electrons density or 
velocity along the radial direction. However, the azimuthal instability in the considered case 
arises due to interaction of the beam with the self-magnetic, rather than external magnetic 
field. 

The azimuthal instability leads to the decrease of the REB critical current due to the 
decrement of the longitudinal velocity, Vz, and the increase of the velocity spread of electrons 
in the region of the formation of the vortex electron structure. As a consequence, VC is 
formed in the region of the vortex structure onset, where the REB space charge density is 
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FIG. 5. Similar to Fig.[3]but for the following parameters: Bq = 30kGs, / = 10 kA; We = 480 keV; 
t2-ti = 0.2 ns 

maximal due to a large number of electrons with low energy. 

Note, generally the vortex structure may be formed in the system without the VC onset. 
The azimuthal instability is developed in the REB when its current is greater than a certain 
threshold value. So, if the REB current exceeds this threshold value and it is less than the 
critical current (when the VC is observed in the system) only the vortex structure is formed 
in the system. 

With the growth of the external magnetic field value the azimuthal instability becomes 
to be suppressed by the focusing force of the external magnetic field which is directed 
oppositely to the centrifugal force. Fig. S] corresponding to the case of Bq ~ -Bmax (^^^ 
Fig. |2^) demonstrates clearly this effect. Actually, the rotating vortex electron structure 
(see Fig. [3f) and d) disappears in the drift space and the electron beam fills the whole space 
without the essential azimuthal inhomogeneity. The space charge density decreases and, 
as a consequence, the critical current for the VC formation increases in this case. So, the 
suppression of the REB azimuthal instability in the applied external magnetic field leads to 
the growth of the REB critical current with the increase of the external magnetic field Bq. 
Note, for the ultra-relativistic electron beams with energy We > 0.7 MeV this increase starts 
simultaneously with the growth of the magnetic field i?o > (curves 3 and 4 in Fig. [2^), 
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whereas for the beams with the less energy {We ~ (0.4 -f- 0.7) MeV) the increase of the REB 
critical current starts for the values of the external magnetic field Bq > Blf^ (see curves 1 
and 2 in Fig. 

A further growth of the magnetic field {Bq > -B^ax) leads to the decrease of the REB 
critical current value and its saturation at constant level for Bq ~ -B-mm similarly to the 
weakly relativistic case^^*^"^. Such behavior of Icr{Bo) is determined by the mechanism of 
the charged particle transversal dynamics limitation by means of the external magnetic field 
(see Fig. E]). As a consequence, the space charge density of the REB also increases with the 
growth of the external magnetic field, and the critical beam current decreases. When the 
value of the external magnetic field induction equals to -B^in REB transversal dynamics 
appears to be completely suppressed. Fig. [5] shows that the transversal beam dynamics is 
not observed in the system in the case of the strong magnetic field Bq > -B^j„. 

The mentioned difference in the curves behavior in Fig. [2] (cf. curves 1 and 2 with curves 3 
and 4) for the case of weak external magnetic fields < Bq < B\ is the consequence 
of the competition of two described above physical mechanisms occurring in the system 
with the growth of the external magnetic field. The first mechanism (the suppression of 
the azimuthal instability of the REB) leads to the critical current value increase. The 
second one, which is connected with the charged particle transversal dynamics limitation 
in external focusing magnetic field, conducts on the contrary to the decrease of critical 
current. For the electron beams with energy We ~ (0.4 0.7) MeV (see curves 1 and 2 in 
Fig. [2^) the both mechanisms have an important role: the second mechanism determines 
the behavior of the dependencies for small values of external magnetic field < Bq < 5*^, 
and the first one begins to dominate for B^j^ < -Bo < -Bmax- the case of ultra-relativistic 
electron beams with energy We > 0.7 MeV when self-magnetic fields are much stronger, 
the formation of VC is determined mainly by the mechanism of vortex electron structure 
formation, whereas the suppression of REB transversal dynamics with the growth of external 
magnetic field has a minor role. So, the first physical mechanism dominates the second one 
considerably. As a consequence the critical current dependencies in this case (see curves 3 
and 4 in Fig. [2^) demonstrate the immediate increase with the external magnetic field growth 
(0 < i?o < Bi^^^). 

Taking into account that the external magnetic field Bq > -B^m is strong enough to 
neglect the influence of the self-magnetic fields, the value -B^in be easily estimated 
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analytically. Let the REB with current Icr has the radius Rb at the injection plane and the 
character radius Rye — at the VC region {Rye > Rb a result of Coulomb's repulsion of 
electrons of beam; the value of Rye is taken from the simulation). Moving in the system 
between the points with radii Rb and Rye in the presence of the external magnetic field, the 
electrons acquire the angular momentum. This momentum is proportional to the difference 
of induction fluxes across the cross-sections of the REB at the points with radii Rb and Rye, 
respectively^^: 

where dO/dt is the azimuthal velocity of electrons. The motion of electrons of REB is 
determined by the action of centrifugal force = jomer [dO /dt)'^ , Coulomb's repulsion 
force Fk = —eE^ and Lorentz force Fl = —er{dd/dt)B (here e and rrie are the charge and 
the mass of the electron, respectively, r is the radial coordinate of electron, Er is the radial 
component of the space charge field intensity). Then, one can write the motion equation 
for the boundary electron of the beam, taking into account the above, the relation and 
equation d^r/dt^ = {2'qVo / 'y q) d'^ r /dz'^: 



d'r , r]B^ 

-riyc 



1 - 



Rb \ 

Rye J 



ATlSo^V^'^RyC 



0, (5) 



dz^ 8\/o7o 

where Vq is the accelerating voltage. 

Eq. ([5]) implies that there is such external magnetic field, -B^m; which REB keeps the 
constant radius in the system. Actually, if we put (P'rjdz^ = in Eq. ^ (it means the lack 
of the acceleration in the radial direction), we obtain the quadratic equation for B^. The 
solution of this equation gives the required value of the external magnetic field -Bmm when 
the dependency /cr(-Bo) saturates (see Fig. [2]): 



-BL.v, — Ryc^ 



mm 



\ireoifl'^V%(Kc-Rt)' 



Fig. [6] demonstrates analytical curve (|6]) (solid line) and the numerically calculated values 
(points) of the external magnetic field -B^m different beam energies We- One can 

see that relation ([6]) describes the results of the numerical simulation correctly and gives 
the accurate values of the magnetic field -B^m being close to the corresponding numerical 
results. The value -B^j„ increases with the growth of the beam energy. We, due to the 
following reason. We would remind that -Bmm is such external magnetic field value for 
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FIG. 6. The comparison of the external magnetic field value B^^^^ defined according to Eq. ([6]) 
(solid curve) with the value calculated with the help of 3D numerical simulation (dots) for 

annular REB 

which the transversal dynamics of electron beam is completely suppressed in the system. 
Obviously, the electron inertia increases with the growth of the beam energy, We, due to 
the electron velocities increase and, as a consequence, the relativistic growth of the electron 
mass is observed. As a result, the greater external magnetic field value is required for the 
electrons in the system to be focused (the limitation of the beam transversal dynamics), 
therefore, Bmin increases monotonously with the growth of the REB energy. The behavior 
of the dependency BminiWe) on the energy We is quahtatively similar to the weakly and 
strong relativistic cases (see Fig. [6]). 

IV. EFFECT OF THE ELECTRON BEAM THICKNESS ON THE 
VIRTUAL CATHODE FORMATION CONDITIONS 

In the previous sections we have analyzed the dynamics of a thin annular REB with a 
fixed thickness, d (see Section [Tll)- Let us consider here the results of numerical simulation 
of the influence of the beam thickness, d, on the critical current of REB as well as on the 
dynamics of the vortex structure in the system. 

Fig. [7] illustrates the dependency of the REB critical current on the beam thickness 
d. One can see that the critical current value decreases monotonously with the growth of 
d. More homogeneous filling of the drift tube by the electron beam with the increasing 
of the beam thickness results in greater potential sagging in the beam drift space. As a 
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FIG. 7. Dependency of the REB critical current on the normalized beam thickness d/Ri, for the 
external magnetic field induction value Bq = 3kGs and the electron beam initial energy value 
We = 480 keV 

consequence, the beam critical current decreases with the growth of the beam thickness that 
agrees qualitatively with the analytical dependency of the critical current value on the beam 
thickness in the case of ID motion of electron beam ^r^. 





FIG. 8. Projections of the instantaneous positions of charged particles of the solid electron beam at 
the transverse plane (r, 0) for Zg = 9.5 mm at the consecutive time points ti and t2 {t2 — ti = 0.2 ns) 
for Bq = 3kGs (upper figures) and Bq = 30kGs (lower figures); I = 7.5 kA, We = 480 keV. Only 
particles behind the projection plane {z > Zg) are shown in Figures 
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Fig. [8] shows the typical projections of the instantaneous positions of charged particles of 
the solid REB at the transverse plane (r, 6) for the beam current being greater the critical 
value and the different characteristic values of the external magnetic field, Bq. Analyzing 
these figures, one can see that the vortex structure in the solid REB is formed (see Fig. [H]a,&) 
similarly to the case of the annular electron beam in the weak external magnetic fields (see 
Section IIIIBj compare Fig. [8]a,6 and Fig. [3]). Let us note that the obtained results for 
the solid REB agree qualitatively with the experimental investigation of the intense REB 
carried out in the Ref.~ where the tracings of witness place damage patterns (see Fig. 6 (left 
panel) in Ref.~) demonstrate the presence of the vortex structure in the solid beam in the 
case of weak external focusing magnetic field. The same result was obtained for different 
beam thicknesses. In particular, with increase of external magnetic field we observe the 
suppression of the vortex pattern formation (see Fig. [8]c,c?) which is also in good agreement 
with an experiment (see Fig. 6 (right panel) in Ref.'^). 

V. CONCLUSION 

Summarizing the obtained results, it is obvious that the VC formation and electron 
structure dynamics in the REB differ considerably in the cases of the weakly and ultra 
relativistic beams. We shown that the nonlinear dynamics of the REB in the weak external 
magnetic fields is defined by the self-magnetic field of REB which leads to the azimuthal 
instability and the vortex electron structure formation. Such behavior reduces the value of 
the REB critical current. The growth of the external magnetic field causes the suppression 
of the azimuthal instability and, consequently, the critical current increases. In the case of 
the strong external magnetic field the critical current value decreases and saturates similarly 
to the weakly relativistic case^"^. The beam thickness does not effect considerably on the 
process of the vortex structure formation, however, the growth of the beam thickness leads 
to the decrease of the beam critical current. 
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